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1. Introduction 

In the last few years, enormous progress has been made in understanding noncritical 
string theory. One line of development started with the work of Refs. fl], |2|, |3| , in the 
context of D-branes of Liouville theory. These and subsequent works were inspired 
by the beautiful CFT computations that gave convincing evidence for the consistency 
of these branes(|, ^ ||, as well as Sen's picture of the decay of unstable D-branes via 
tachyon condensation^]. Another independent line of development that has proved 
important was the attempt to formulate new matrix models to describe noncritical 
string theories and their deformations, including black hole deformations^], [| [IUL O]. 

Some of the important new results are related to nonperturbatively stable type 
fermionic strings ]T^|, |T3 , but even in the bosonic context, many old and new puzzles 



concerning matrix models as well as Liouville theory have been resolved. For c < 1 
matter coupled to Liouville theory, a beautiful picture emerged of a Riemann surface 
governing the semiclassical dynamics of the model. Both ZZ and FZZT branes were 
identified as properties of this surface: the former are located at singularities while 



the latter arise as line integrals. This picture was obtained in Ref. fl4}| within the 
continuum Liouville approach and subsequently re-derived in the matrix model for- 
malism in Ref. [13] using earlier results of Ref. [T^J. However, later it was realised| 
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that the exact, as opposed to semiclassical, picture is considerably simpler: the Rie- 
mann surface disappears as a result of Stokes' phenomenon and is replaced by a 
single sheet. In the exact (quantum) case, correlation functions of macroscopic loop 
operators go from multiple-valued functions to the Baker- Akhiezer functions of the 
KP hierarchy, which are analytic functions of the boundary cosmological constant. 
Thus, for these models (and also their type extensions) a rather complete picture 
now exists. 

Another remarkable development in this context is an explicit proposal to un- 
derstand open-closed string duality starting from open string field theory. This was 
presented in Ref. ||18|| and implemented there for the (2, q) series of minimal models 
coupled to gravity (which can be thought of as perturbations of the "topological 
point" or (2,1) minimal model). The basic idea of Ref. ||18|| was to evaluate open 
string field theory on a collection of N FZZT branes in the (2, 1) closed string back- 
ground. This leads to the Kontsevich matrix model |Tj|, which depends on a constant 
matrix A whose eigenvalues are the N independent boundary cosmological constants 
for this collection of branes. Now the Kontsevich model computes the correlators 
of closed- string observables in the same (2, 1) background. So this relationship was 
interpreted as open-closed duality, following earlier ideas of Sen[pU|. 

A different way of understanding what appears to be the same open-closed 
duality emerged in Ref. ||17|| for the (2,1) case. Extending some older observa- 



tions in Ref. ||21||, the authors showed that if one inserts macroscopic loop operators 
det(xi — $), representing FZZT branes (each with its own boundary cosmological 
constant Xi) in the Gaussian matrix model, and takes a double-scaling limit, one 
obtains the Kontsevich matrix model. The constant matrix A in this model again 
arises as the boundary cosmological constants of the FZZT branes 1 . 

The situation is more complicated and less well-understood for c = 1 matter 
coupled to Liouville theory, namely the c = 1 string. The results of FZZT were 
derived for generic Liouville central charge c^, but become singular as cl — > 25, the 
limit that should give the c = 1 string. Attempts to understand FZZT branes at c = 1 
(Refs. |25|, 26| |) rely on this limit from the c < 1 case which brings in divergences and 
can therefore be problematic. In particular, there is as yet no definite computation 
exhibiting open-closed duality at c = 1 starting from open string field theory in the 
c = 1 Liouville background. One should expect such a computation to give rise to 
the c = 1 analogue of the Kontsevich matrix model, namely the Kontsevich-Penner 
model 2 of Ref. 0. 

In the present work we take a different approach to understand D-branes and 
open-closed duality in the c = 1 string, more closely tied to the approach of Refs. [117, 



1 This has been generalised |Z| by starting with macroscopic loops in the double-scaled 2-matrix 
models that describe (p, 1) minimal model strings. After double-scaling, one obtains the generalised 
Kontsevich models of Refs.[^3|, 

2 This model is valid only at the selfdual radius R = 1. 
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22|j . The obvious point of departure at c = 1 would be to consider macroscopic 
loops in the Matrix Quantum Mechanics (MQM) and take a double-scaling limit. 
Indeed, FZZT branes at c = 1 have been investigated from this point of view, for 



example in Refs. |28|, 29]. However, we will take an alternative route that makes use 
of the existence of the Normal Matrix Model (NMM) |1TJ for c = 1 string theory (in 
principle, at arbitrary radius R). This model is dual in a certain precise sense to 
the more familiar MQM, namely, the grand canonical partition function of MQM is 
the partition function of NMM in the large- N limit. Geometrically, the two theories 
correspond to different real sections of a single complex curve. More details about 



the interrelationship between MQM and NMM can be found in Ref.[ll 



One good reason to start from the NMM is that it is a simpler model than MQM 
and does not require a double-scaling limit. Also, it has been a longstanding question 
whether the KP model and NMM are equivalent, given their structural similarities, 
and if so, what is the precise map between them. It is tempting to believe that 
open-closed duality underlies their mutual relationship. Indeed, the NMM does not 
have a parameter suggestive of a set of boundary cosmological constants, while the 
KP model has a Kontsevich-type constant matrix A. So another natural question is 
whether the eigenvalues of A are boundary cosmological constants for a set of FZZT 
branes/macroscopic loop operators of NMM. 

In what follows we examine these questions and obtain the following results. First 
of all we find a precise map from the NMM (with arbitrary tachyon perturbations) 
to the KP model, thereby demonstrating their equivalence. While the former model 
depends on a non-Hermitian matrix Z constrained to obey [Z, Z^\ = 1, the latter 
is defined in terms of a positive definite Hermitian matrix M. We find that the 
eigenvalues Zj and are related by m, = z^. The role of the large- N limit in the 
two models is slightly different: in the KP model not only the random matrix but 
also the number of parameters (closed string couplings) is reduced at finite N. On 
the contrary, in the NMM the number of parameters is always infinite for any N, but 
one is required to take N —>■ oo to obtain the right theory (this was called "Model 
I" in Ref . ||1 If ) . The two models are therefore equivalent only on a subspace of the 



parameter space at finite N, with the limit iV — > oo being required to obtain full 
equivalence. This is an important point to which we will return. 

Next in §[| we consider macroscopic loop operators of the form det(£ — Z) in the 
NMM, and show that these operators when inserted into the NMM, decrease the value 
of the closed-string tachyon couplings in a precise way dictated by the Kontsevich- 
Miwa transform. On the contrary, operators of the form 1/ det(£ — Z) play the role of 
increasing, or turning on, the closed-string tachyon couplings. In particular, insertion 
of these inverse determinant operators in the (partially unperturbed) NMM leads to 
the Kontsevich-Penner model. (By partially unperturbed, we mean the couplings 
of the positive-momentum tachyons are switched off, while those of the negative- 
momentum tachyons are turned on at arbitrary values.) Calculationally, this result 
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is a corollary of our derivation of the KP model from the perturbed NMM in §|j. 

These results bear a rather strong analogy to the emergence of the Kontsevich 
model from the insertion of determinant operators at c < 1[|17|. In both cases, 



the parameters of macroscopic loop operators turn into eigenvalues of a Kontsevich 
matrix. Recall that in Ref . |17| , one inserts n determinant operators into the N x N 
Gaussian matrix model and then integrates out the Gaussian matrix. Taking iV — > oo 
(as a double-scaling limit) we are then left with the Kontsevich model of rank n. In 
the c = 1 case, we insert n inverse determinant operators in the NMM. As we will 
see, N — n of the normal matrix eigenvalues then decouple, and we are left with a 
Kontsevich-Penner model of rank n (here one does not have to take iV — > oo). We 
see that the two cases are rather closely analogous. 



The main difference between our case at c = 1 and the c < 1 case of Ref. 17 



that we work with inverse determinant rather than determinant operators. However 
at infinite n we can remove even this difference: it is possible to replace the inverse 
determinant by the determinant of a different matrix, defining a natural pair of 
mutually "dual" Kontsevich matrices 3 . In terms of the dual matrix, one then recovers 
a relation between correlators of determinants (rather than inverse determinants) and 
the KP model. 

In the concluding section we examine a peculiar property of the NMM, namely 
that it describes the c = 1 string even at finite N, if we set N = u, where v is the 



analytically continued cosmological constant v = —ifi. This was noted in Ref. [[11 



where this variant of the NMM was called "Model II" . Now it was already observed 
in Ref.f2~?| that setting N = u in the KP model (and giving a nonzero value to one of 
the deformation parameters) reduces the KP model to the original Kontsevich model 
that describes (2, q) minimal strings. Thus we have a (two-step) process leading from 
the NMM to the Kontsevich model. However, we also know from Ref.[|n| that the 
Kontsevich model arises from insertion of macroscopic loops in the double-scaled 
Gaussian matrix model. We will attempt to examine to what extent these two facts 
are related. 



2. Normal Matrix Model 

We start by describing the Normal Matrix Model (NMM) of c = 1 string theory HIT | 



and making a number of observations about it. The model originates from some well- 
known considerations in the Matrix Quantum Mechanics (MQM) description of the 
Euclidean c = 1 string at radius R. Here, R = 1 is the selfdual radius, to which we 
will specialise later. The MQM theory has discrete "tachyons" 7^, of momentum 
where k e Z. Let us divide this set into "positive tachyons" T^, k > and "negative 



3 This dual pair is apparently unrelated to the dual pair of boundary cosmological constants at 
c < 1. 
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tachyons" Xfc, k < 0. (The zero-momentum tachyon is the cosmological operator and 
is treated separately). We now perturb the MQM by these tachyons, using coupling 
constants t k ,k > for the positive tachyons and tk, k > for the negative ones. 

The grand canonical partition function of MQM is denoted Z(/j,,tk,tk)- At t k = 
tk = 0, it can easily be shown to be: 

z(im, t k = o, t k = o) = J] r ( — ^1 - z/i + - J (2.i) 

n>0 ^ ' 

But this is also the partitition function of the matrix integral: 
Znmm = [ [dZdZ^} e-« w V>*) 

[dZdZ^] e tK-K^ t ) H +[|(fl-i)+(^-^)] logzzt) (2 ^ 2) 

where z/ = — z/i and iV — > oo. Here Z, Z^ are N x N matrices satisfying: 

[Z, Z ] ] = (2.3) 



Since the matrix Z commutes with its adjoint, the model defined by Eq. ( |2.2j ) is 
called the Normal Matrix Model (NMM) 4 . 

The equality above says that the unperturbed MQM and NMM theories are 
equivalent. The final step is to note that the tachyon perturbations correspond to 
infinitely many Toda "times" in the MQM partition function, which becomes a r- 
function of the Toda integrable hierarchy. The same perturbations on the NMM side 
are obtained by adding to the matrix action the terms: 

oo 

W(Z, Zt) - W(Z, Z^) + uJ2 (tkZ k + hZ ]k ) (2.4) 

fc=i 

It follows that the Normal Matrix Model, even after perturbations, is equivalent to 
MQM. 

The equivalence of the full perturbed MQM and NMM gives an interesting inter- 
pretation of the perturbations in NMM in terms of the Fermi surface of the MQM. 
The unperturbed MQM Hamiltonian is given by: 

Hq = -ti(-h 2 ^- - X 2 ) (2.5) 
2 v dX 2 ' K ' 



4 For the most part we follow the conventions of Ref.^lJ. However we use the transcription 
(l/ifi.)them — * Vms, and /ithcm — * lus- The partition function depends on the ratio (/J,/ih)them — * 
v us = —i/j, us . Our conventions for the NMM will be seen to match with the conventions of Ref.|27| 
for the KP model. Note that the integral is well-defined for all complex v with a sufficiently large 
real part, ft can then be extended by analytic continuation to all complex values of the parameter 
z/, other than those for which the argument of the T function is a negative integer. This is sufficient, 
since everything is ultimately evaluated at purely imaginary values of v. 
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where X is an N x N Hermitian matrix (here the compactification radius is R). 
In the 5'L r (A^)-singlet sector this system is described by N non-relativistic fermions 
moving in an inverted harmonic oscillator potential. The eigenvalues of X describe 
the positions of these fermions. In terms of eigenvalues the Hamiltonian can be 
written as: 

1 N 

ffo-E^-^ ( 2 - 6 ) 

Pi being the momenta conjugate to X{. We now want to consider perturbations of 
Eq. ( |2.6| ) by tachyon operators. For this it is convenient to change variables from p, 
x to the "light cone" variables x±: 

x ± = (2-7) 

Since [p,x] = —ih it follows that = —ih also. The MQM Hamiltonian in 

terms of the new variables is: 

ihN 



H = - ^ x+iX-i — (2.i 



i=l 

In the phase space (x + ,x_) the equation of the Fermi surface for the unperturbed 
MQM is given by: 

= fi (2.9) 

The tachyon perturbations to the MQM Hamiltonian Hq are given in terms of the 
new variables by: 

N 

k>l i=l 

In the above equation the v's are determined in terms of the fs from the orthonor- 
mality of the Fermion wavef unctions. The conventions chosen above simplifies the 
connection with NMM perturbations. The Fermi surface of the perturbed MQM is 
given by: 

x + x- = fi + ^2 (kt± k x£ + v± k x± R j (2.11) 

k>l 

The equivalence between NMM and MQM relates the tachyon perturbations 
in Eq. ( ^.4[ ) and Eq. ( |2.1(J| ) with the following identification between the tachyon 
operators of the two models: 

n 

trX* = trZ n 
trX% = trZ tn 
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The coefficients t± are the same as t, t in the NMM. This means that any tachyon 
perturbation in the NMM is mapped directly to a deformation of the Fermi surface 
of MQM by Eq. (|2TTT1 ). 

At the self dual radius R=l, the NMM simplifies and the full perturbed partition 
function can be written as: 

z NMM (t,t) = J [dzdtf] ^(-^Hv-mio^-^ (fczMizt*)) (2 12) 

We note several properties of this model. 

(i) The unperturbed part depends only on the combination ZZ^ and not on Z, Z^ 
separately. 

(ii) The model can be reduced to eigenvalues, leading to the partition function: 

/N 
i=l 

(2.13) 

(iii) The model is symmetric under the interchange tk «-> tk, as can be seen by 
interchanging Z and Z\ In spacetime language this symmetry amounts to the trans- 
formation X — > —X where X is the Euclidean time coordinate, which interchanges 
positive and negative momentum tachyons. 

(iv) The correlator: 

(trZ fel trZ fc2 ■ • • tiZ km tiZ ]l \iZ^ 2 ■ ■ ■ tiZ^ n ) tk= - tk=0 (2.14) 

vanishes unless 

m n 

This correlator is computed in the unperturbed theory. The above result follows by 
performing the transformation: 

Z -> e i0 Z (2.16) 

for some arbitrary angle 9. The unperturbed theory is invariant under this trans- 
formation, therefore correlators that are not invariant must vanish. In spacetime 
language this amounts to the fact that tachyon momentum in the X direction is 
conserved. 

(v) As a corollary, we see that if we set all t k = 0, the partition function becomes 
independent of %: 

Znmm{®-> tk) — Znmm(Q, 0) (2-17) 

(vi) For computing correlators of a finite number of tachyons, it is enough to 
turn on a finite number of £fc, tk, i.e. we can always assume for such purposes that 
tk,tk = 0, k > k max for some finite integer k max . In that case, apart from the log 
term we have a polynomial matrix model. 
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(vii) We can tune away the log term by choosing v — N. This choice has 



been called Model II in Ref.|TT|. In this case the model reduces to a Gaussian 
model (but of a normal, rather than Hermitian, matrix) with perturbations that are 
holomorphic + antiholomorphic in the matrix Z (i.e., in the eigenvalues Zj). If we 
assume that the couplings tk, tk vanish for k > k max , as in the previous comment, then 
the perturbations are also polynomial. We will return to this case in a subsequent 
section. 



3. The Kontsevich-Penner or Wqq model 

The Kontsevich-Penner or model |27j (for a more detailed review, see Ref.|30 



is a model of a single positive-definite hermitian matrix, whose partition function is 
given by: 

(A,t) = {detA) v I [dM] e H-» MA +(»- N ^ M -»^=^ Mk ) (3.1) 



<£>KP\ 

where tk are the couplings to negative-momentum tachyons, N is the dimensionality 
of the matrix M and A is a constant matrix. The eigenvalues of this matrix determine 
the couplings £& to positive-momentum tachyons via the Kontsevich-Miwa (KM) 
transform: ^ 

t k = — r tT ( A ' k ) ( 3 - 2 ) 

UK 



This model is derived by integrating the equations found in Ref.[|3T]. The 
parameter v appearing in the action above is related to the cosmological constant \x 
of the string theory by v = —ip,. The model can also be obtained from the Penner 
matrix model[|3^, [33[ after making a suitable change of variables (as explained in 
detail in Ref.|3(J) and adding perturbations. 

We now note some properties that are analogous to those of the NMM, as well 
as others that are quite different. 

(i) By redefining MA — > M we can rewrite the partition function without any 
factor in front, as: 

Z KP (A,t) = J [dM] e ^M + (^N)io S u-,E^UMA-^) (3 3) 

(ii) This model has no radius deformation, and describes the c = 1 string theory 
directly at self dual radius R = 1. 

(iii) In view of the logarithmic term, the model is well-defined only if the integral 
over the eigenvalues m; of the matrix M is restricted to the region > 0. 

(iv) The model can be reduced to eigenvalues, leading to the partition function: 



Z KP (A,t) = (f[a)j I f[dmt ^^eE^iC-Wi+^iog^ES^W) (3.4) 
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(vi) In the representation Eq. (|3.1|) , the operators trM fc describe the negative- 
momentum tachyons. But there are no simple operators that directly correspond to 
positive-momentum tachyons. Nevertheless this model generates tachyon correlators 
of the c = 1 string as follows: 

(T kl T k2 ■ ■■T km T^ h ■ . .T.O = ■ ■ ■ -i^lk^ " ' W n hgZKP (3 ' 5) 



where derivatives in t k are computed using Eq. (3^) and the chain rule. 

(v) The symmetry of the partition function under the interchange of t k ,t k is not 
manifest, since one set of parameters is encoded through the matrix A while the 
other appears explicitly. 

(vi) The transformation 

A — > aA, t k ^a k t k (3.6) 

for arbitrary a, is a symmetry of the model (most obvious in the representation 
Eq. ( |3.3|) ). As a consequence, the tachyon correlators satisfy momentum conserva- 
tion. 

(vi) The partition function satisfies the "puncture equation" : 

Z KP (A -e,t k + 5 kA e) = e evH ^Z KP (A,t k ) (3.7) 
as can immediately be seen from Eq. 



4. Equivalence of the matrix models 

4.1 N = 1 case 

We start by choosing the self dual radius R = 1, and will later comment on what 
happens at other values of R. As we have seen, in the perturbed NMM there are two 
(infinite) sets of parameters t k} t k , all of which can be chosen independently. This is 
the case even at finite N, though the model describes c = 1 string theory only at 
infinite N (or at the special value N = u, as noted in Ref.|TT|], a point to which we 



will return later). In contrast, the Kontsevich-Penner model has one infinite set of 
parameters t k , as well as N additional parameters from the eigenvalues of the matrix 
A. The latter encode the t k , as seen from Eq. ( |3.2|) above. From this it is clear that 
at finite N, there can only be independent parameters t k (k = 1, 2, . . . N) while 
the remaining ones (t k , k > N) are dependent on these. 

This makes the possible equivalence of the two models somewhat subtle. To 
understand the situation better, let us compare both models in the limit that is 
farthest away from N —>■ oo, namely N — 1. While this is a "toy" example, we will 
see that it provides some useful lessons. 
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In this case the NMM partition function is: 

Z NM M, N =i(t k ,tk) = f dzdz e --+(-i)io g --E£i(*^ fe +^^) (4.1) 

while the Kontsevich-Penner partition function is: 

Z KPiN=1 (a,t k ) =a u J dm e -^+(^)^-^T=i~t^ k ( 4 . 2 ) 

We will now show that the two integrals above are equivalent if we assume that t k 
in the NMM is given by: 

t k = -\a- k (4.3) 

UK 



which is the KM transform Eq. ( |3.2| ) in the special case where A is a 1 x 1 matrix, 
denoted by the single real number a. Note that this determines all the infinitely 
many t k in terms of a. 

To obtain the equivalence, insert the above relation and also perform the change 
of integration variable: 

z = ^/me id (4.4) 
in the NMM integral. Then we find that (up to a numerical constant): 

Z NM M,N=x(a,t k ) = I dmdd e — +(- 1 ) lo s-+S£ 1 i(^) fe ^-E^ 1 4(v^)'= e -- 

dmdO l—- e - vm lo s™-^^i4(v^) fc e-^ ( 45 ) 

-/me ta 



1 - 



Strictly speaking the last step is only valid for -s/m/a < 1, since otherwise the infinite 
sum fails to converge. Hence we fix m and a to satisfy this requirement and continue 
by evaluating the ^-integral. This can be evaluated by defining e~ ld = w and treating 
it as a contour integral in w. We have 

d9 %=r- s - -> dw 1 —= (4.6) 

1 — W — - — 

a a 

Since the rest of the integrand is well-defined and analytic near w = 0, we capture the 
simple pole at w — y/m/a. That brings the integrand to the desired form. Now we 
can lift the restriction ^fmja < 1, and treat the result as valid for all m by analytic 
continuation. Therefore we find: 

^NMM,N=i{a,tk) — / am e 6 ^ k - 1 

= Z KPN=1 {a,t k ) (4.7) 
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Thus we have shown that the perturbed lxl Normal Matrix Model at R — 1 is 
equivalent to the lxl Kontsevich-Penner model. However, this equivalence only 
holds when we perform the lxl KM transform, which fixes all the perturbations t& 
in terms of a single independent parameter a (while the tk are left arbitrary). 

An important point to note here is the sign chosen in Eq. ((Of). Changing the sign 
(independently of k) amounts to the transformation tk — > —tk- This is apparently 
harmless, leading to some sign changes in the correlation functions, but there is no 
way at N = 1 (or more generally at any finite N) to change a (or the corresponding 
matrix A) to compensate for this transformation. The sign we have chosen, given 
the signs in the original NMM action, is therefore the only one that gives the KP 
model. This point will become important later on. 

Returning to the NMM-KP equivalence at N — 1, it is interesting to generalise 
it by starting with the NMM at an arbitrary radius R instead of R = 1 as was the 
case above. As seen from Eq. Q2.2p , the coupling of the log term is modified in this 
case as: 

(i/-l)-±(J*-l) + (J&/-l) (4.8) 

and also the bilinear term zz is modified to (zz) R . The above derivation goes through 
with only minor changes, and we end up with: 

Z N MM,N=i(a,t k ) = a***- 1 *" J dm e -K»n«)"+[i(Ji-i)+(^-i)] logm-, ESLi?^ (49) 

This appears to suggest that there is a variant of the Kontsevich-Penner model valid 
at arbitrary radius (or at least arbitrary integer radius, since otherwise it may become 
hard to define the integral). This would be somewhat surprising as such a model 
has not been found in the past. As we will see in the following subsection, the above 
result holds only for the N = 1 case. Once we go to N x N matrices, we will see 
that NMM leads to a KP matrix model only at R = 1, consistent with expectations. 

Another generalisation of the above equivalence seems more interesting. In prin- 
ciple, even for the lxl matrix model, we can carry out a KM transform using an 
n x n matrix A where n is an arbitrary integer. Indeed, there is no logical reason 
why the dimension of the constant matrix A must be the same as that of the ran- 
dom matrices occurring in the integral. The most general example of this is to take 
N x N random matrices Z, in the NMM and then carry out a KM transform with 
A being an n x n matrix. The "usual" transform then emerges as the special case 
n = N . Of course all this makes sense only within the NMM and not in the KP 
model. If n ^ N then the KP model, which has a trMA term in its action, cannot 
even be defined. So we should not expect to find the KP model starting with the 
NMM unless n = N, but it is still interesting to see what we will find. 

Here we will see what happens if we take N = 1 and n > 1. The full story 
will appear in a later subsection. Clearly the KM transform Eq. ( |3.2j ) permits more 
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independent parameters tk as n gets larger. Let us take the eigenvalues of A to be 
a±, 02, ... , a n . Then it is easy to see that: 

7 (n 7 \- I rim rift p -™+(^-l)log™+ELiE£i iff l^'^^E^i^tV^)^-* 9 

j nr=ia-^) 

Converting to the w variable as before, we now encounter n poles. Picking up the 
residues, we get: 



J NMM,N- 



4oi3k) = ( dm e -^ + (,-i)io g m^| 1 

J l=1 aj 



1 -"E2Li?*(^)* 



(4.11) 

This in turn can be expressed as a sum over n 1 x 1 Kontsevich-Penner models: 

" 1 

Znmm,n=i(pu tk) = Y\ 7i ~^ZKP,N=i{au tk) (4-12) 



1=1 n^/ 

Note that if in this expression we take a n —>■ oo, one of the terms in the above equa- 
tion (corresponding to I = n) decouples, and a n also drops out from the remaining 
terms. Therefore we recover the same equation with n — ► n — 1. In this way we can 
successively decouple all but one of the aj's. 

To summarise, at the level of the lxl NMM, we have learned some interesting 
things: this model is equivalent to the lxl KP model if we specialise the parameters 
tk to a 1-parameter family via the KM transform, while it is equivalent to a sum over 
n different lxl KP models if we specialise the parameters tk to an n-parameter 
family. We also saw a 1 x 1 KP model arise when we are at a finite radius R ^ 1. 
In the next section we will see to what extent these lessons hold once we work with 
N x N random matrices. 

4.2 General case 

In this section we return to the N x N Normal Matrix Model. With the substitution 
Eq. ( |3.2| ) (where A is also an N x N matrix), its partition function becomes: 

z NMM = J {dZdtf)A -^ t +(-^)^^ t +2:£ 1 i^)^-Er =1 *^t'=) (4 13) 

or, in terms of eigenvalues: 

N 

i~Zi + {v-N) J2fLl lo S z iZi 



Znmm = f f[d 2 Zl A(z)A(z)e-»^^ 



x e w=1 * =1 *W e -"E£iE2Li**3f (4.14) 
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where A(^) is the Vandermonde determinant. Because of the normality constraint 
[Z, Z*] = there is only one Vandermonde for Zi and one for z\. 

The sum over k in the second line of Eq. (|4.14 ) converges if ^ < 1 for all i,j, in 
which case it can be evaluated immediately giving: 



N N 



Z NMM = f f[d 2 Zl \A(z)\ 2 f[ r 4 1I e S£i[-^+(«-^^-"5:£= 1 i*2f] (4.15) 



1 = 1 jj = l 



To make contact with the Penner model, first change variables z% — > y/mit 
and then replace e~ l9i by Wi as before. Then we get d 2 Zi — >■ drrii — and: 



N „ N , N 

~'NMM 



AT 

x TT 1 e Ef=i[-^rn i +(u-N)\ogm i -u^ =1 t k ( y /nHW i ) k ] /^gx 

The contour integrals can be evaluated once this is rewritten in the more convenient 
form: 

N „ V TV 



Z N mm = / J^drrii ® J^rfwj (y/nTiWj - v /m~w i ) {^/mlwi - ^/rfTjWj) 

J i=l ^ i=l i<j 

x JT^ c ^, J Ii [-^m^+(^-V)logm.i--t/^gl 1 t fc (Vrn~Wi) fc ] (4 17) 



Next we pick up the residues at the poles. During the intermediate steps, we 
will assume that the eigenvalues aj of the matrix A are non-degenerate. From the 
above expression, each integration variable Wi has a pole at each of the values: 



Wi = ^ (4.18) 



for all j. Thus the contributions can be classified by the set of poles: 

(wi, w 2 ,..., w N ) = I , , . . . , J (4.19) 

\ aj 1 Oj 2 aj N j 

We now notice that the set (ji, J2, • • • , Jjv) must consist of distinct elements, in other 
words it forms a permutation of (1,2, . . . ,N). This is because if two values of ji 
coincide, one of the Vandermonde factors of the type (y/mlwj — ^rrijWi) vanishes 
and there is no contribution. 

We start by considering the simplest permutation, the identity, namely: 

tii,h,...,j N ) = (l,2,...,N) (4.20) 
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In this case the residues from the denominator and Vandermonde factors become: 



/ ^/mimj ^/m~m~\ /rrii mA -q 1 1 1^. ,('"'"/ _ m j a i 



. dj di j V a i a i J f — v ^ A(a 

while the exponential measure factor becomes: 



e 



-jmn+^-JV) log mi-vj:T=i **> (^) ' 



(4.21) 



(4.22) 



It is easy to check that for all the other possible permutations of (ji, j 2 , . . . , Jn) 
besides the identity permutation, a corresponding permutation of the integration 
variables 777.4 brings the above answer (exponential measure as well as prefactors) 
back to the same form as for the identity permutation. This means that (dropping 
a factor of j^) we have proved: 

N yjN 



Ui<j( m i a j -rrijai) £^ J-^+^-jv) lo gmi -, ££° =1 (^) * 
A(a) 6 



Znmm = / ] J <^ m i 
i=i 

= (tl^Y / 11^ ^eSf^h^^+^^^-E^^] (4.23) 

4=1 ^ 4=1 V / 

where in the last step we have replaced m« — > mjOj. 

This is precisely the eigenvalue representation Eq. ( |3.4| ) of the KP matrix model 
Eq. Thus we have provided a direct proof of equivalence of the perturbed 

Normal Matrix Model and the Kontsevich-Penner model. Notice that in performing 
the KM transform we reduced the independent t k of the NMM to a finite number, 
namely N, so that eventually the N — > oo limit is required in order to encode all the 
independent parameters. 

In the previous subsection we considered taking different ranks for the constant 
matrix A arising in the KM transform and the random matrix Z. The most general 
case is to take Z to be N x N and A to be n x n. The computation is a simple 
extension of the one done above. We find the following results. When n > N we 
again get a sum over Kontsevich-Penner models. The number of terms in the sum is 
the binomial coefficient n C^. This is a generalisation of the result given in Eq. ( |4.12j ) 
for N — 1, where we found n terms. In the general case let us denote by a,{i t q the 
z th element of the set formed by one possible choice of N a^s from a total of n, the 
index I labeling the particular choice. The complementary set, formed by the rest of 
the cii's is denoted by n, the index i taking n — N values. We then have: 

n C N N N-n ^ 

z NM M{ahh) = y~] rr rr -? — - — \Z KP (ain,t k ) (4.24) 

1=1 4=1 - =1 II— —— 

1 1 V a {i,i} 
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so that the NMM is again expressed as a sum over KP models. 

The other case, n < N, can be obtained by starting with n = N and successively 
decoupling N — n eigenvalues Oj by taking them to infinity. This is similar to what 
we observed in the N = 1 case following Eq. ( 4.12| ). In the present case one can 



easily show that N — n matrix eigenvalues rrii also decouple in this limit (apart from 
a normalisation). In fact, it is straightforward to derive the formula: 

lim Z^\A<"Kt k ) = r( ^ + 1) Z&-^-'U) (4-25) 

which can then be iterated. Thus after N — n eigenvalues cij are decoupled, we find up 
to normalisation the KP model of rank n. As we remarked in the introduction, this 
exhibits a strong analogy to the insertion of n determinant operators in the Gaussian 
model, as described in Ref. Ill 711, where the result is the n x n Kontsevich model 5 . 



4.3 Radius dependence 

Finally, we can ask what happens to the radius-dependent NMM under the above 
procedure. Again the steps are quite straightforward and one arrives at the following 



generalisation of Eq. ( 4.23 ): 



J NMM,R 



V \{R-\)+v f" , Aim) 



A(a) 

x e T,ti[-Hmiai) R +[^(R-l)+^R-N)] Iogmi-i/E2Li **"»?] (4.26) 

The problem is that the above eigenvalue model cannot (as far as we can see) be 
converted back to a matrix model. The key to doing so in the R = 1 case was the 
linear term Yli m i a i m the action, which (after absorbing the Vandermondes and 
using the inverse of the famous Harish Chandra formula) can be summed back into 
tvMA. The quantity Y2i( m i a i) R cannot be converted back into a matrix trace unless 
R = l. 

This clarifies a longstanding puzzle: while a KP model could only be found at 
R = 1, the NMM exists and describes the c = 1 string for any R. We see now that 
the correct extension of the KP model to R ^ 1 is the eigenvalue model given by 
Eq. ( |4.26|) above, but unfortunately this does not correspond to a matrix model. 



5. Loop operators in the NMM 

In this section we will examine loop operators in the NMM. Our goal here is to 
understand whether correlation functions of these operators can be related to the 

5 This paragraph corrects an error in a previous version of this paper. As a result, the analogy 
with c < 1 is now stronger than we had previously claimed. We are grateful to the referee for 
helpful suggestions in this regard. 
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Kontsevich-Penner model of Ref.[27], thereby providing the c = 1 analogue of the 
corresponding observations in Refs. |17|, Though there are some similarities, we 
will also find some striking differences between this and the c < 1 case. 

Macroscopic loops in a model of random matrices $ are described by insertions 
of the operator: 

W{x) = tr log(x - $) (5.1) 

which creates a boundary in the world sheet. Here x is the boundary cosmo logical 
constant. The corresponding generating function for multiple boundaries is|3^, [35, 
3§,§7|,[T7|: 

e w & = det(x - $) (5.2) 

Such operators have been studied extensively in c < 1 matrix models, describing 
(p, q) minimal models coupled to 2d gravity. 

We will consider expectation values of operators of the form det(a — Z) in the 
NMM, where a is a real parameter. These operators create a hole in the dual graph 
in the Feynman diagram expansion of the matrix model. Since the NMM has vertices 
that are holomorphic/antiholomorphic in Z, the dual graph will have faces that are 
dual to Z or Z^ . The loop operator det(a — Z) creates a hole in a Z-face, while its 
complex conjugate creates a hole in a Z^-face. 

As we would expect, this means that the correlators are complex, but we have 
the identity 6 : 

(jTdet(oi-Z)) _ = (Y[det(a t -&))_ (5.3) 

\ / t k ,t k \ / tk,tk 

i i 

where on the RHS the role of the deformations tk, tk has been interchanged. Therefore 
as long as we consider correlators only of det(aj — Z) or det(aj — Z') the result is 
effectively the same. As we will see in a moment, a stronger statement is true: on the 
subspace of parameter space dictated by the KM transform, the unmixed correlators 
are individually real. Later we will also consider mixed correlators. 
As a start, notice that in the lxl case, 



Z NMM ,N=i{t k = 0,t k ) = I d 2 z e -»**H»-D****-» 

d 2 z (a - z)—^— e~ vzi H-("-i)ks**-" ET=J^ k 
(a — z) 

d 2 z (a - z) e -^^M^)^gzz-uY, k < L 1 (t k z k +t k z k ) 

~((a-z)) Z NMMN=1 (t° k ,t k ) (5.4) 
a\ /t° k ,t k 



6 Here and in the rest of this section, all correlators are understood to be normalised correlators 
in the NMM. 
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where the expectation value in the last line is evaluated in the NMM with 

t° k = -\a- k (5.5) 

UK 

We see that the t° k dependence drops out in the RHS because insertion of the loop 
operator cancels the dependence in the partition function. In fact, more is true: even 
the t k dependence cancels out between the different factors on the RHS. This is a 
consequence of the property exhibited in Eq. ( |2.17|) . 
A more general statement in the lxl case is: 

2>NMM,N=i(tk — t% tk) — —((a — z)\ _ Z NMMtN= i(t k ,tk) (5.6) 

In other words, insertion of the macroscopic loop operator has the effect of decreasing 
the value of tk, leaving t k unchanged. 

In the more general case of N x N random matrices, the corresponding result is 
as follows. The expectation value of a single exponentiated loop operator det (a — Z) 
is: 

det(a-Z)) _ = — a (5.7) 

I t k ,t k ■£>NMM{tk,tk) 

with t° k again given by Eq. (|5.5| ). Now we would like to consider multiple loop 
operators. Therefore consider the expectation value: 

n 

JJdet(fli-Z)\ _ (5.8) 



As noted in Ref. ||17||, this can be thought of as a single determinant in a larger space. 
Define the nxn matrix A = diag(ai, a,2, ■ ■ ■ , a n ) and extend it to an (n + N) x (n + N) 
matrix A ® IjvxAt- Similarly, extend the N x N matrix Z to an (n + N) x (n + N) 



matrix l nxn ® Z. Now we can write 



n n N 

I J det(ai -Z) = det(A ®1-1®Z) = JJ {^[(a, - z s ) (5.9) 

i=l i=l j=l 

Rewriting this as: 

n 

Y[ det(oi -Z) = (det A) N det(l ® 1 - A" 1 ® Z) (5.10) 



and expanding the second factor, we find: 

i=1 *k.*fc _ Z NMM {t kl t 
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where now: 

4 = -±-trA~ k (5.12) 

Thus we see that macroscopic loop correlators in this model are obtained by simply 
shifting the parameters tk in the partition function, the shift being given by the KM 
transform. 

The above considerations can be extended to mixed correlators as follows. Con- 
sider correlation functions of the form: 



JJdet(ai-Z) jjdet(6 i - Z ] )\ (5.13) 



Then, defining the m x m matrix B = diag(6i, 62, • • • , b m ), the parameters t° as in 



Eq. (|5.12| ), and the parameters by: 

t = -\trB- k (5.14) 
vk 

we find 

TT det(a, - Z) TT det(6, - #)) = (det A det b) n 

i=1 . =1 ' t k ,t k ^NMM\tk-,tk) 

(5.15) 

In the above we have seen how to re-express correlations of loop operators in 
terms of shifted closed-string parameters. This in itself is quite reminiscent of an 
open-closed duality. However we did not yet encounter the KP model. To do so, we 
note that besides the exponentiated loop operator det (a — Z), we can consider its 
inverse: 1/ det (a — Z). Just as insertion of det (a — Z) has the effect of decreasing 
each tk by t° k given by Eq. (|5.5| ), insertion of the inverse operator increases by the 
same amount. 

Thus we may consider correlators like: 

lidettfli-Z)) = (det A)" (det(l®ll-A- 1 ®Z) N ) ^' 16 ^ 



i=i 



As before, the two factors of the direct product in the above equation refer tonxii 
and N x N matrices. It is easy to see that the correlation function on the RHS has 
the effect of increasing the tk by t° k as given in Eq. ( |5.12|) . 



Although in principle n and iV are independent, here we will consider the case 
n = N . Now the inverse operator 

' ' (5-17) 



det(l®]l-A- 1 ®Z) 
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has already made an appearance in §|4], where one finds it in the eigenvalue basis (see 
for example Eq. ( |4.15|) ): 

N 1 

n ( 5 - ig ) 

i,j=l a j 



The interesting property of the inverse determinant operators is that they can 
be used to create the KP model starting from the partially unperturbed NMM (where 
tk = but tk are arbitrary). Computationally this is similar to the derivation in §^ 
of the KP model from the perturbed NMM. Thus we have: 

Z NMM (0,tk) = Z KP (A,tk) (5.19) 

o,i fc 



det(l ® 1 - A- 1 ® Z) 



Here 2 NMM (0,tk) can be replaced by Znmm(0, 0) as we have noted previously. This 
equation then is the precise statement of one of our main observations, that inverse 
determinant expectation values in the (partially unperturbed) NMM give rise to the 
KP partition function. 

It is clearly desirable to have a target space interpretation for these loop oper- 
ators. Since the NMM is derived from correlators computed from matrix quantum 
mechanics, in principle one should be able to understand the loop operators of NMM 
starting from loop operators (or some other operators) in MQM. While that is beyond 
the scope of the present work, we will instead exhibit some suggestive properties of 
our loop operators and leave their precise interpretation for future work. 

In matrix models for the c < 1 string, which are described by constant random 
matrices, exponentiated loop operators are determinants just like the ones discussed 
here for the NMM. In those models it has been argued that the loop operators rep- 
resent FZZT branes. One striking observation is that in the Kontsevich/generalised 
Kontsevich description of c < 1 strings, the eigenvalues of the constant matrix A 
come from the boundary cosmological constants appearing in the loop operators. 
Moreover, Eq. (|5.2| ) has been interpreted as evidence that the FZZT-ZZ open strings 
there are fermionic P?|, |TT |. 



In the present case, we see that the parameters a* in the loop operators turn 
precisely into the eigenvalues of the constant matrix A of the Kontsevich- Penner 
model. We take this as evidence that our loop operators are likewise related in some 
way to FZZT branes. Indeed, one is tempted to call them FZZT branes of the NMM. 
Pursuing this analogy further, the role played by inverse determinants in the present 
discussion appears to suggest that the corresponding strings in the NMM are bosonic 
rather than fermionic. But the relationship of these operators to the "true" FZZT 
branes of matrix quantum mechanics remains to be understood, as we have noted 
above 7 . 



7 In light of the discussions about the MQM Fermi surface in we can give an interpretation 
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In the limit of infinite N, the inverse loop operators depending on a matrix A can 
be thought of as loop operators for a different matrix A. Thus, only in this limit, the 
inverse determinant operators can be replaced by more conventional determinants. 
This proceeds as follows. We have already seen that the KM transformation Eq. ( |3.2j ) 
encodes infinitely many parameters tk via a constant N x N matrix A, in the limit 
N —* oo. Now for fixed tk, suppose we considered the (very similar) transform: 

t k = ^-tTA- k (5.20) 

UK 

that differs only by a change of sign. The point is that this apparently harmless 
reversal of the tk brings about a significant change in the matrix A. Moreover this 
is possible only in the infinite N limit, since we are trying to satisfy: 

trA~ k = -trA- k (5.21) 

for all k. Now it is easy to see that the matrices A and A satisfy the following 
identity: 

det(l ® H - A' 1 ® Z) = (5.22) 

det(l ® 1 - A- 1 ® Z) 

Therefore a correlator of inverse loop operators can be rewritten in terms of usual 
loop operators using: 

j * " —(det(A®H-l®Z)) (5.23) 



det(A® 1 - 1®Z)/ (detAA) 

In the light of our previous observation that inverse determinant operators might 
indicate the bosonic nature of FZZT-ZZ strings at c = 1, it is tempting to think of 
Eq. (|5.23|) as a statement of fermi-bose equivalence! 

In terms of the operator det(A (g) 1 — 1® Z), we can make the statement that its 
insertion into the partially unperturbed NMM gives rise to the KP model depending 
on the "dual" Kontsevich matrix A. 



6. Normal matrix model at finite N 

The correspondence between NMM and KP model demonstrated in §|] is valid for 
any N, as long as the parameters of the former are restricted to a subspace. The 
NMM itself is supposed to work at iV — > oo, in which case this restriction goes 



away. However, as noted in Ref . ||1 1|| , there is another way to implement the NMM: 



to both determinant and inverse determinant operators in the NMM. Since their insertions lead to 
opposite shifts in the tk's, by virtue of the equivalence between MQM and NMM discussed above 
we can map each one directly to a corresponding deformation to the Fermi surface, which can be 



read off from Eq. (2.11). This fact should facilitate direct comparison with the MQM. 
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by setting iV = vR (which amounts to iV = v for R — 1), which they labelled as 
"Model II". In other words, these authors argue that: 



lim Z NMM (N, t, v) = Z NMM (N = vR, t, v) (6.1) 

Thus the NMM describes the c = 1 theory at this finite value of N, after analytically 
continuing the cosmological constant fi = iv to an imaginary value 8 . 

The key property of this choice is that the logarithmic term in the matrix po- 
tential of the NMM gets tuned away. Let us take R = 1 from now on. Suppose 
we evaluate the expectation value of the inverse determinant operator at this N (for 
the moment we assume that this special value is integral). For N insertions of the 
inverse determinant, it gives the KP model with N = v. Thus, as one would expect, 
the log term of the KP model is also tuned away. Now if we choose tk = c5ks, with 
c some constant, then the KP model reduces to the Kontsevich model, as observed 
in Ref . [f2?H . This shows that the Kontsevich model is a particular deformation of the 
c = 1 string theory after analytic continuation to imaginary cosmological constant 
and condensation of a particular tachyon (T 3 ) . Note that at the end of this procedure, 
the rank of the Kontsevich matrix is the same as that of the NMM matrix. 

As mentioned earlier, there is a different route to the Kontsevich model starting 
from the Gaussian Matrix Model (GMM) Jl7[]. Here one starts with a Gaussian matrix 
model of rank N, with iV insertions of the determinant operator, and takes N — > oo 
as a double- scaling limit by focussing on the edge of the eigenvalue distribution. The 
result is the Kontsevich matrix model. This time the rank N of the original matrix 
has disappeared from the picture (it was sent to infinity) while the Kontsevich matrix 
inherits its rank from the number of determinant insertions N. 

A diagram of the situation is given in Fig.l. ^From the figure one sees that 
the diagram can be closed if we find a suitable relation of the NMM to the Gaussian 
matrix model. This is not hard to find at a qualitative level. In fact with v = N 
and tk,tk = the NMM is a Gaussian matrix model. We choose the rank to be 
N. The NMM eigenvalue distribution p(z, z) is constant inside a disc in the z-plane 



(for R = l)fllT|. If we look at a contour along the real axis in the z-plane, then the 



effective eigenvalue distribution 

p(x) = I dy p{x,y) (6.2) 



is a semi-circle law, and we find the Gaussian matrix model. However, this picture of 
eigenvalue distributions is valid only at large N. Inserting N determinant operators 
and taking N — > oo as a double-scaling limit, one recovers the Kontsevich model. In 



8 Whereas the authors of Ref. 1 11 presented this as the analytic continuation of N to the imaginary 
value — ijji, we prefer to think of it as continuing the cosmological constant [i to the imaginary value 
iN. 
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NMM 



N x N 



N 1 \ 
j2i det(a,-Z)'0,t k 



KP 



N 

( n detfaj - M)) 
i=i 



N xN 
v = N 

tk = <5fe,3 



GMM - Kontsevich 

N -> oo 

N x N N x N 

Figure 1: Two routes from NMM to Kontsevich 



this way of proceeding, the cubic coupling of the Kontsevich model is switched on 
automatically during the double-scaling limit. In the alternative route through the 
KP model, one has to switch on the coupling t 3 by hand. A more detailed under- 
standing of these two routes and their relationship should illuminate the question of 
how minimal model strings are embedded in c = 1. We leave this for future work. 



7. Conclusions 

We have established the equivalence between two matrix models of the c = 1 string 
(at selfdual radius): the Normal Matrix Model of Ref.JTTJ and the Kontsevich- Penner 



model of Ref. ||27f. Both matrix models were initially found as solutions of a Toda 
hierarchy, so this equivalence is not very surprising. However, it is still helpful to 
have an explicit derivation, which also uncovered a few subtleties. Also we ended up 
showing why the KP matrix model does not exist at radius R ^ 1. 

The more interesting aspect of this equivalence is that correlation functions of 
inverse determinant operators in the partially unperturbed NMM give rise to the 
KP model. This is analogous to corresponding results in Refs. [TT|, |2"2"fl, with two im- 



portant differences. In those cases, one considered determinants rather than inverse 
determinants, and their correlators were computed in a double-scaled matrix model. 
In the NMM there is no double-scaling as it already describes the grand canonical 
partition function of the double-scaled Matrix Quantum Mechanics. Another differ- 
ence is that the N of the final (KP) model is equal to that of the NMM, and part 
of the matrix variables in NMM survive as the matrices of the KP model. All this 
suggests that, if one makes an analogy with the topological minimal models, the 
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NMM occupies a position half-way between the original matrix model arising from 
dynamical triangulation of random surfaces (which requires a double-scaling limit 
to describe continuum surfaces) and the final "topological" model. If this is true, 
we may have only described half the story of open-closed duality at c = 1 while 
the correspondence between MQM and NMM constitutes the previous half. Further 
work may lead to a more coherent picture of the steps involved and thereby a deeper 
understanding of open-closed string duality at c = 1. 

As we commented earlier, the inverse determinant operators seem to suggest 
bosonic statistics for FZZT-ZZ branes (at least in the NMM context) in contrast to 
fermionic statistics for c < 1. Another way to think of this is that both determinant 
and inverse determinant operators expand out to give the same set of macroscopic 
loops, the only difference being a minus sign for an odd number of loops in the latter 
case. Alternatively one can think of the basic loop operator as being changed by a 
sign to — trlog(a — Z). Either way, the role of inverse determinant operators clearly 
calls for further investigation. 

We commented earlier that trying to take the c = 1 limit of c < 1 FZZT 
correlators is problematic and therefore a derivation of the KP model from open- 
string field theory analogous to Ref. [IS] has not been forthcoming. While this may 
yet be achieved, the situation recalls a historical parallel. In the 1990's, attempts 
to derive c = 1 closed string theory as a limit of the c < 1 theories were not very 
successful. Eventually it was found that at least at selfdual radius, the c = 1 string is 
a nonstandard case - rather than a limit - of the c < 1 models. This was understood 
by going over to the topological |38|, |39], ^(J rather than conventional, formulation of 
these string theories. It emerged that while the (p, q) minimal models for varying 
q were described by topological models labelled by an integer k = p — 2 > (for 
example, SU(2))./U(1) twisted Kazama-Suzuki models or twisted Af = 2 Landau- 
Ginzburg theories with superpotential X k+2 ), the c = 1 string at selfdual radius 



was instead described by "continuations" of these models to k — —3 41, 40, 42, 43 



rather than the more naive guess one might have made, namely k — > oo. Therefore 
progress on FZZT branes at c = 1 in the continuum formulation might most naturally 
emerge in the context of topological D-branes in the twisted SU{2)-$/U{l) Kazama- 
Suzuki model or X^ 1 Landau-Ginzburg theory. Indeed, Ref.|0| represents important 
progress in this direction, and the Kontsevich model has been obtained there in the 
topological setup, predating the more recent derivations of Refs. ||TE|, [17). In fact, the 
KP model of Ref. |27|| was also obtained in Ref.[44|. 

Extension of the NMM/KP models to include winding modes of the c = 1 string, 



as well as a better understanding of 2d black holes from matrix models || |9|, [10], |TT 
remain open problems and perhaps the open-closed duality studied here will be 
helpful in this regard. 

We have not pursued here an observation made in Ref. [p0[1 that the KP model 
simplifies when we exponentiate the matrix variable via M = e*. The resulting 



- 23 - 



model, which was named the "Liouville matrix model" there, is suggestive of iV 
D-instantons moving in a Liouville plus linear potential. A similar exponentiation 
can be carried out in the NMM. In either case this is an almost trivial change of 
variables, therefore it does not seem important for the considerations in the present 
paper. However, in the light of the present work, these changes of variables might 
lead to new and more satisfying interpretations of the matrix models themselves. 

As a final comment, we note that open-closed duality has in recent times been 
given a more fundamental basis in the Gopakumar programme [|45|, |46| , 47|j where 
closed string theory is proposed to be derived from quite general large- N field the- 
ories. Now this programme is expected to apply not just to noncritical strings but 
to all string theories. We know that the Kontsevich and Penner models compute 
topological invariants of the moduli space of Riemann surfaces, but the above works 
seem to suggest that these models play a role in more complicated string theories 
too. If so, equivalences and open-closed dualities such as we have discussed here may 
have more far-reaching implications than just providing examples in simplified string 
backgrounds. 
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